21 results
Stall onset on aerofoils at low to moderately high Reynolds number flows – ADDENDUM
- Wallace J. Morris II, Zvi Rusak
-
- Journal:
- Journal of Fluid Mechanics / Volume 943 / 25 July 2022
- Published online by Cambridge University Press:
- 21 June 2022, E1
-
- Article
- Export citation
Transonic flows of single-phase supercritical fluids over thin airfoils
- Zvi Rusak, Akashdeep Singh Virk
-
- Journal:
- Journal of Fluid Mechanics / Volume 915 / 25 May 2021
- Published online by Cambridge University Press:
- 17 March 2021, A61
-
- Article
- Export citation
-
A small-disturbance model for a steady, two-dimensional, inviscid and transonic flow of a single-phase real gas around a thin airfoil is presented. The approach explores the nonlinear interactions among near-sonic speed of the flow, small thickness ratio of the airfoil and upstream properties of the fluid. The gas thermodynamic properties are related by a general equation of state. Information about thermodynamic modelling of the gas is lumped into one similarity parameter, $K_G$, related to the fundamental derivative of gas dynamics. The flow field is described by a modified transonic small-disturbance problem. The theory applies to any working fluid of interest. Model problems are derived for steam flows described by the perfect, van der Waals, virial and Redlich–Kwong gas equations of state. Predictions are compared according to the various gas models under various free-stream operating conditions from low subcritical to high supercritical thermodynamic states to gain insights into the sensitivity of the small-disturbance problem solution to thermodynamic modelling of the gas. Results show that transonic flows are independent of gas modelling at low subcritical thermodynamic conditions. However, at near-critical and supercritical thermodynamic conditions, transonic flow behaviour is significantly sensitive to gas modelling and variations of $K_G$. The upstream flow critical Mach number increases as the flow approaches thermodynamic critical state and a wider range of upstream Mach numbers can be found where pressure drag is zero. However, at supercritical conditions, $K_G$ increases, resulting in lower critical Mach numbers and higher pressure drags.
Vortex flows of moist air with non-equilibrium and homogeneous condensation
- Zvi Rusak, Gerald A. Rawcliffe, Yuxin Zhang
-
- Journal:
- Journal of Fluid Mechanics / Volume 885 / 25 February 2020
- Published online by Cambridge University Press:
- 07 January 2020, A38
-
- Article
- Export citation
-
A small-disturbance model is presented for the complex dynamics of vortex flows of moist air in a straight, circular pipe with non-equilibrium and homogeneous condensation. The model explores the nonlinear interactions among the vortex near-critical swirl ratio and the small amount of water vapour in the air. The condensation rate is calculated according to classical nucleation and droplet growth models. The asymptotic analysis gives the similarity parameters that govern the flow problem. These are the flow inlet swirl ratio $\unicode[STIX]{x1D714}$, the inlet Mach number $Ma_{0}$, the initial humidity $\tilde{\unicode[STIX]{x1D714}}_{0}$, the number of water molecules in a characteristic fluid element $n_{C}$, the inlet centreline super-saturation ratio $S_{0}$ and the ratio of characteristic condensation and flow time scales $K$. Also, the flow field may be described by an ordinary first-order nonlinear differential equation for the flow evolution coupled with a set of four first-order ordinary differential equations along the pipe for the calculation of the condensate mass fraction. An iterative numerical scheme which combines the Runge–Kutta integration technique for the flow dynamics with Simpson’s integration rule for the calculation of the condensation variables is developed. Specifically, equilibrium states are determined, including the possibility of the appearance of multiple states under the same boundary conditions, and the stability characteristics of these states are described. The model is used to study the effects of humidity and of energy supply from nanoscale condensation processes on the large-scale dynamics of vortex flows as well as the effect of flow swirl on condensation processes in swirling flows.
Near-sonic pure steam flow with real-gas effects and non-equilibrium and homogeneous condensation around thin airfoils
- Akashdeep Singh Virk, Zvi Rusak
-
- Journal:
- Journal of Fluid Mechanics / Volume 884 / 10 February 2020
- Published online by Cambridge University Press:
- 13 December 2019, A30
-
- Article
- Export citation
-
A small-disturbance asymptotic model is derived to describe the complex nature of a pure water vapour flow with non-equilibrium and homogeneous condensation around a thin airfoil operating at transonic speed and small angle of attack. The van der Waals equation of state provides real-gas relationships among the thermodynamic properties of water vapour. Classical nucleation and droplet growth theory is used to model the condensation process. The similarity parameters which determine the flow and condensation physics are identified. The flow may be described by a nonlinear and non-homogeneous partial differential equation coupled with a set of four ordinary differential equations to model the condensation process. The model problem is used to study the effects of independent variation of the upstream flow and thermodynamic conditions, airfoil geometry and angle of attack on the pressure and condensate mass fraction distributions along the airfoil surfaces and the consequent effect on the wave drag and lift coefficients. Increasing the upstream temperature at fixed values of upstream supersaturation ratio and Mach number results in increased condensation and higher wave drag coefficient. Increasing the upstream supersaturation ratio at fixed values of upstream temperature and Mach number also results in increased condensation and the wave drag coefficient increases nonlinearly. In addition, the effects of varying airfoil geometry with a fixed thickness ratio and chord on flow properties and condensation region are studied. The computed results confirm the similarity law of Zierep & Lin (Forsch. Ing. Wes. A, vol. 33 (6), 1967, pp. 169–172), relating the onset condensation Mach number to upstream stagnation relative humidity, when an effective specific heat ratio is used. The small-disturbance model is a useful tool to analyse the physics of high-speed condensing steam flows around airfoils operating at high pressures and temperatures.
Swirling flow states of compressible single-phase supercritical fluids in a rotating finite-length straight circular pipe
- Nguyen Ly, Zvi Rusak, Shixiao Wang
-
- Journal:
- Journal of Fluid Mechanics / Volume 849 / 25 August 2018
- Published online by Cambridge University Press:
- 21 June 2018, pp. 576-614
-
- Article
- Export citation
-
Steady states of inviscid, compressible and axisymmetric swirling flows of a single-phase, inert, thermodynamically supercritical fluid in a rotating, finite-length, straight, long circular pipe are studied. The fluid thermodynamic behaviour is modelled by the van der Waals equation of state. A nonlinear partial differential equation for the solution of the flow streamfunction is derived from the fluid equations of motion in terms of the inlet flow specific total enthalpy, specific entropy and circulation functions. This equation reflects the complicated, nonlinear thermo-physical interactions in the flows, specifically when the inlet state temperature and density profiles vary around the critical thermodynamic point, flow compressibility is significant and the inlet swirl ratio is high. Several types of solutions of the resulting nonlinear ordinary differential equation for the axially independent case describe the flow outlet state when the pipe is sufficiently long. The approach is applied to an inlet flow described by a solid-body rotation with uniform profiles of the axial velocity and temperature. The solutions are used to form the bifurcation diagrams of steady compressible flows of real fluids as the inlet swirl level and the centreline inlet density are increased at a fixed inlet Mach number and temperature. Focus is on heavy-molecule fluids with low values of $R/C_{v}$. Computed results provide theoretical predictions of the critical swirl levels for the exchange of stability of the columnar state and for the appearance of non-columnar states and of vortex breakdown states as a function of inlet centreline density. The difference in the dynamical behaviour between that of a calorically perfect gas and of a real gas is explored. The analysis sheds new fundamental light on the complex dynamics of high-Reynolds-number, compressible, subsonic swirling flows of real gases.
Dynamics of a perturbed solid-body rotation flow in a finite-length straight rotating pipe
- Chunjuan Feng, Feng Liu, Zvi Rusak, Shixiao Wang
-
- Journal:
- Journal of Fluid Mechanics / Volume 846 / 10 July 2018
- Published online by Cambridge University Press:
- 16 May 2018, pp. 1114-1152
-
- Article
- Export citation
-
Direct numerical simulations are used to study the three-dimensional, incompressible and viscous flow dynamics of a base solid-body rotation flow with a uniform axial velocity entering a rotating, finite-length, straight circular pipe. Steady in time profiles of the axial, radial and circumferential velocities are prescribed along the pipe inlet. The convective boundary conditions for each velocity flux component is set at the pipe outlet. The simulation results describe the neutral stability line in response to either axisymmetric or three-dimensional perturbations in a diagram of Reynolds number ( $Re$ , based on inlet axial velocity and pipe radius) versus the incoming flow swirl ratio ( $\unicode[STIX]{x1D714}$ ). This line is in good agreement with the neutral stability line recently predicted by the linear stability theory of Wang et al. (J. Fluid Mech., vol. 797, 2016, pp. 284–321). The computed time history of the velocity components at a certain point in the flow is used to describe three-dimensional phase portraits of the flow global dynamics and its long-term behaviour. They show three types of flow evolution scenarios. First, the Wang & Rusak (Phys. Fluids, vol. 8 (4), 1996, pp. 1007–1016) axisymmetric instability mechanism and evolution to a stable axisymmetric breakdown state is recovered at certain operational conditions in terms of $Re$ and $\unicode[STIX]{x1D714}$ . However, at other operational conditions with same $\unicode[STIX]{x1D714}$ but with a higher $Re$ , a second scenario is found. The axisymmetric breakdown state continues to evolve and a spiral instability mode appears on it and grows to a rotating spiral breakdown state. Moreover, at higher levels of $\unicode[STIX]{x1D714}$ a third scenario is found where there exists a dominant three-dimensional spiral type of instability mode that agrees with the linear stability theory of Wang et al. (J. Fluid Mech., vol. 797, 2016, pp. 284–321). The growth of this mode leads directly to a spiral type of flow roll-up and nonlinearly saturates on a rotating spiral type of vortex breakdown. The Reynolds–Orr equation is used to reveal the mechanism that drives all the instabilities as well as the nonlinear global flow evolution. At high swirl ratios, the confined kinetic energy in the swirling flow can be triggered to be released through various physical agents, such as the asymmetric inlet–outlet conditions, that eliminate axial homogeneity along the pipe and induce flow instabilities and evolution to breakdown states. It is also shown that local instability analysis or its extension using the assumption of a weakly non-parallel flow to conduct convective instability–absolute instability analyses is definitely not related to any of the instability modes found in the present study. Moreover, a stability study based on the strongly non-parallel flow character, including axial inhomogeneity due to a finite-domain boundary conditions, must be conducted to reveal instabilities in such flows.
Swirling flow states in finite-length diverging or contracting circular pipes
- Zvi Rusak, Yuxin Zhang, Harry Lee, Shixiao Wang
-
- Journal:
- Journal of Fluid Mechanics / Volume 819 / 25 May 2017
- Published online by Cambridge University Press:
- 27 April 2017, pp. 678-712
-
- Article
- Export citation
-
The dynamics of inviscid-limit, incompressible and axisymmetric swirling flows in finite-length, diverging or contracting, long circular pipes is studied through global analysis techniques and numerical simulations. The inlet flow is described by the profiles of the circumferential and axial velocity together with a fixed azimuthal vorticity while the outlet flow is characterized by a state with zero radial velocity. A mathematical model that is based on the Squire–Long equation (SLE) is formulated to identify steady-state solutions of the problem with special conditions to describe states with separation zones. The problem is then reduced to the columnar (axially-independent) SLE, with centreline and wall conditions for the solution of the outlet flow streamfunction. The solution of the columnar SLE problem gives rise to the existence of four types of solutions. The SLE problem is then solved numerically using a special procedure to capture states with vortex-breakdown or wall-separation zones. Numerical simulations based on the unsteady vorticity circulation equations are also conducted and show correlation between time-asymptotic states and steady states according to the SLE and the columnar SLE problems. The simulations also shed light on the stability of the various steady states. The uniqueness of steady-state solutions in a certain range of swirl is proven analytically and demonstrated numerically. The computed results provide the bifurcation diagrams of steady states in terms of the incoming swirl ratio and size of pipe divergence or contraction. Critical swirls for the first appearance of the various types of states are identified. The results show that pipe divergence promotes the appearance of vortex-breakdown states at lower levels of the incoming swirl while pipe contraction delays the appearance of vortex breakdown to higher levels of swirl and promotes the formation of wall-separation states.
Near-critical swirling flow of a viscoelastic fluid in a circular pipe
- Zvi Rusak, Nguyen Ly, John A. Tichy, Shixiao Wang
-
- Journal:
- Journal of Fluid Mechanics / Volume 814 / 10 March 2017
- Published online by Cambridge University Press:
- 06 February 2017, pp. 325-360
-
- Article
- Export citation
-
The interaction between flow inertia and elasticity in high-Reynolds-number, axisymmetric and near-critical swirling flows of an incompressible and viscoelastic fluid in an open finite-length straight circular pipe is studied at the limit of low elasticity. The stresses of the viscoelastic fluid are described by the generalized Giesekus constitutive model. This model helps to focus the analysis on low fluid elastic effects with shear thinning of the viscosity. The application of the Giesekus model to columnar streamwise vortices is first investigated. Then, a nonlinear small-disturbance analysis is developed from the governing equations of motion. It reveals the complicated interactions between flow inertia, swirl and fluid rheology. An effective Reynolds number that links between steady states of swirling flows of a viscoelastic fluid and those of a Newtonian fluid is revealed. The effects of the fluid viscosity, relaxation time, retardation time and mobility parameter on the flow development in the pipe and on the critical swirl for the appearance of vortex breakdown are explored. It is found that in vortex flows with either an axial jet or an axial wake profile, increasing the shear thinning by decreasing the ratio of the viscoelastic characteristic times from one (with fixed values of the Weissenberg number and the mobility parameter) increases the critical swirl ratio for breakdown. Increasing the fluid elasticity by increasing the Weissenberg number from zero (with a fixed ratio of the viscoelastic characteristic times and a fixed value of the mobility parameter) or increasing the fluid mobility parameter from zero (with fixed values of the Weissenberg number and the ratio of viscoelastic times) causes a similar effect. The results may explain the trend of changes in the appearance of breakdown zones as a function of swirl level that were observed in the experiments by Stokes et al. (J. Fluid Mech., vol. 429, 2001, pp. 67–115), where Boger fluids were used. This work extends for the first time the theory of vortex breakdown to include effects of non-Newtonian fluids.
On the three-dimensional stability of a solid-body rotation flow in a finite-length rotating pipe
- Shixiao Wang, Zvi Rusak, Rui Gong, Feng Liu
-
- Journal:
- Journal of Fluid Mechanics / Volume 797 / 25 June 2016
- Published online by Cambridge University Press:
- 18 May 2016, pp. 284-321
-
- Article
- Export citation
-
The three-dimensional, inviscid and viscous flow instability modes that appear on a solid-body rotation flow in a finite-length straight, circular pipe are analysed. This study is a direct extension of the Wang & Rusak (Phys. Fluids, vol. 8 (4), 1996a, pp. 1007–1016) analysis of axisymmetric instabilities on inviscid swirling flows in a pipe. The linear stability equations are the same as those derived by Kelvin (Phil. Mag., vol. 10, 1880, pp. 155–168). However, we study a general mode of perturbation that satisfies the inlet, outlet and wall conditions of a flow in a finite-length pipe with a fixed in time and in space vortex generator ahead of it. This mode is different from the classical normal mode of perturbations. The eigenvalue problem for the growth rate and the shape of the perturbations for any azimuthal wavenumber $m$ consists of a linear system of partial differential equations in terms of the axial and radial coordinates ($x,r$). The stability problem is solved numerically for all azimuthal wavenumbers $m$. The computed growth rates and the related shapes of the various perturbation modes that appear in sequence as a function of the base flow swirl ratio (${\it\omega}$) and pipe length ($L$) are presented. In the inviscid flow case, the $m=1$ modes are the first to become unstable as the swirl ratio is increased and dominate the perturbation’s growth in a certain range of swirl levels. The $m=1$ instability modes compete with the axisymmetric ($m=0$) instability modes as the swirl ratio is further increased. In the viscous flow case, the viscous damping effects reduce the modes’ growth rates. The neutral stability line is presented in a Reynolds number ($Re$) versus swirl ratio (${\it\omega}$) diagram and can be used to predict the first appearance of axisymmetric or spiral instabilities as a function of $Re$ and $L$. We use the Reynolds–Orr equation to analyse the various production terms of the perturbation’s kinetic energy and establish the elimination of the flow axial homogeneity at high swirl levels as the underlying physical mechanism that leads to flow exchange of stability and to the appearance of both spiral and axisymmetric instabilities. The viscous effects in the bulk have only a passive influence on the modes’ shapes and growth rates. These effects decrease with the increase of $Re$. We show that the inviscid flow stability results are the inviscid-limit stability results of high-$Re$ rotating flows.
Extension to nonlinear stability theory of the circular Couette flow
- Pun Wong Yau, Shixiao Wang, Zvi Rusak
-
- Journal:
- Journal of Fluid Mechanics / Volume 795 / 25 May 2016
- Published online by Cambridge University Press:
- 19 April 2016, pp. 455-493
-
- Article
- Export citation
-
A nonlinear stability analysis of the viscous circular Couette flow to axisymmetric finite-amplitude perturbations under axial periodic boundary conditions is developed. The analysis is based on investigating the properties of a reduced Arnol’d energy-Casimir function $\mathscr{A}_{rd}$ of Wang (Phys. Fluids, vol. 2, 2009, 084104). A weighted kinetic energy of the perturbation, which has a form of ${\rm\Delta}\mathscr{A}_{rd}$, the difference between the reduced Arnol’d function and its base flow value, is used as a Lyapunov function. We show that all the inviscid flow effects as well as all the viscous-dependent terms that are related to the flow boundaries vanish. The evolution of ${\rm\Delta}\mathscr{A}_{rd}$ depends only on the viscous effects of the perturbation’s dynamics inside the flow domain. The requirement for the temporal decay of ${\rm\Delta}\mathscr{A}_{rd}$ leads to two novel sufficient conditions for the nonlinear stability of the circular Couette flow in response to axisymmetric perturbations. The linearized version of these conditions for infinitesimally small perturbations recovers the recent linear stability results by Kloosterziel (J. Fluid Mech., vol. 652, 2010, pp. 171–193). By examining the nonlinear stability conditions, we establish a definite operational region of the viscous circular Couette flow that is independent of the fluid viscosity. In this region of operation, the flow is nonlinearly stable in response to perturbations of any size, provided that the initial total circulation function is above a minimum level determined by the operational conditions of the base flow. Comparisons with historical studies show that our results shed light on the experimental measurements of Wendt (Ing.-Arch., vol. 4, 1933, pp. 577–595) and extend the classical nonlinear stability results of Serrin (Arch. Rat. Mech. Anal., vol. 3, 1959, pp. 1–13) and Joseph & Hung (Arch. Rat. Mech. Anal., vol. 44, 1971, pp. 1–22). When the flow is nonlinearly stable and evolves axisymmetrically for all time, then it always decays asymptotically in time to the circular Couette flow determined uniquely by the set-up of the rotating cylinders. Finally, we derive upper-bound estimates on the decay rate of finite-amplitude perturbations for the solid-body rotation flow between two coaxial rotating cylinders and for the circular Couette flow. We demonstrate via numerical simulations that the theoretical upper bound is relevant to the dynamics of various axisymmetric perturbations tested, where it is strictly obeyed. This present study provides new physical insights into a classical flow problem that was studied for many decades.
Vortex breakdown in premixed reacting flows with swirl in a finite-length circular open pipe
- Zvi Rusak, Jung J. Choi, Nicholas Bourquard, Shixiao Wang
-
- Journal:
- Journal of Fluid Mechanics / Volume 793 / 25 April 2016
- Published online by Cambridge University Press:
- 22 March 2016, pp. 749-776
-
- Article
- Export citation
-
A global analysis of steady states of low Mach number inviscid premixed reacting swirling flows in a straight circular finite-length open pipe is developed. We focus on modelling the basic interaction between the swirl and heat release of the reaction. For analytic simplicity, a one-step first-order Arrhenious reaction kinetics is considered in the limit of high activation energy and infinite Peclet number. Assuming a complete reaction with chemical equilibrium upstream and downstream of the reaction zone, a nonlinear partial differential equation is derived for the solution of the flow stream function downstream of the reaction zone in terms of the specific total enthalpy, specific entropy and circulation functions prescribed at the inlet. Several types of solutions of the nonlinear ordinary differential equation for the columnar flow case describe the outlet states of the flow in a long pipe. These solutions are used to form the bifurcation diagram of steady reacting flows with swirl as the inlet swirl level is increased at a fixed heat release from the reaction. The approach is applied to two profiles of inlet flows, the solid-body rotation and the Lamb–Oseen vortex, both with constant profiles of the axial velocity, temperature and mixture reactant mass fraction. The computed results provide theoretical predictions of the critical inlet swirl levels for the appearance of vortex breakdown states and for the size of the breakdown zone as a function of the inlet flow swirl level, Mach number and heat release of the reaction. For the inlet solid-body rotation, flow is decelerated to breakdown as the inlet swirl is increased above the critical swirl level, and there is a delay in the appearance of breakdown with the increase of the heat release of the reaction. For the inlet Lamb–Oseen vortex at low values of heat release, the critical swirl for breakdown is decreased with the increase of heat release while, at high values of heat release, the appearance of breakdown is delayed to higher incoming flow swirl levels with the increase of heat release. The analysis sheds light on the global dynamics of low Mach number reacting flows with swirl and vortex breakdown and on the interaction between vortex breakdown and heat release that affects the shape of the reaction zone in the domain.
Vortex breakdown of compressible subsonic swirling flows in a finite-length straight circular pipe
- Zvi Rusak, Jung J. Choi, Nicholas Bourquard, Shixiao Wang
-
- Journal:
- Journal of Fluid Mechanics / Volume 781 / 25 October 2015
- Published online by Cambridge University Press:
- 16 September 2015, pp. 3-27
-
- Article
- Export citation
-
A global analysis of steady states of inviscid compressible subsonic swirling flows in a finite-length straight circular pipe is developed. A nonlinear partial differential equation for the solution of the flow stream function is derived in terms of the inlet flow specific total enthalpy, specific entropy and circulation functions. The equation reflects the complicated thermo–physical interactions in the flows. Several types of solutions of the resulting nonlinear ordinary differential equation for the columnar case together with a flow force condition describe the outlet state of the flow in the pipe. These solutions are used to form the bifurcation diagram of steady compressible flows with swirl as the inlet swirl level is increased at a fixed inlet Mach number. The approach is applied to two profiles of inlet flows, solid-body rotation and the Lamb–Oseen vortex, both with a uniform axial velocity and temperature. The computed results provide for each inlet flow profile theoretical predictions of the critical swirl levels for the appearance of vortex breakdown states as a function of the inlet Mach number, suggesting that the results are robust for a variety of inlet swirling flows. The analysis sheds light on the dynamics of compressible flows with swirl and vortex breakdown, and shows the delay in the appearance of breakdown with increase of the inlet axial flow Mach number in the subsonic range of operation. The present theory is limited to axisymmetric dynamics of swirling flows in pipes where the wall boundary layer is thin and attached and does not interact with the flow in the bulk.
An active feedback flow control theory of the axisymmetric vortex breakdown process
- Zvi Rusak, Joshua Granata, Shixiao Wang
-
- Journal:
- Journal of Fluid Mechanics / Volume 774 / 10 July 2015
- Published online by Cambridge University Press:
- 15 June 2015, pp. 488-528
-
- Article
- Export citation
-
An active feedback flow control theory of the axisymmetric vortex breakdown process in incompressible swirling flows in a finite-length straight circular pipe is developed. Flow injection distributed along the pipe wall is used as the controller. The flow is subjected to non-periodic inlet and outlet conditions where the inlet profiles of the axial velocity, circumferential velocity and azimuthal vorticity are prescribed, along with no radial velocity at the outlet. A long-wave asymptotic analysis at near-critical swirl ratios, which involves a rescaling of the axial distance and time, results in a model problem for the dynamics and the nonlinear control of both inviscid and high-Reynolds-number ($\mathit{Re}$) flows. The approach provides the bifurcation diagram of steady states and the stability characteristics of these states. In addition, an energy analysis of the controlled flow dynamics suggests a feedback control law that relates the flow injection to the evolving maximum radial velocity at the inlet. Computed examples of the flow dynamics based on the full Euler and Navier–Stokes formulations at various swirl levels demonstrate the evolution to near-steady breakdown states when swirl is above a critical level that depends on $\mathit{Re}$. Moreover, applying the proposed feedback control law during flow evolution shows for the first time the successful and robust elimination of the breakdown states and flow stabilization on an almost columnar state for a wide range of swirl (up to at least 30 %) above critical. The feedback control cuts the natural feed-forward mechanism of the breakdown process. Specifically, in the case of high-$\mathit{Re}$ flows, the control approach establishes a branch of columnar states for all swirl levels studied, where in the natural flow dynamics no such states exist. The present theory is limited to the control of axisymmetric flows in pipes where the wall boundary layer is thin and attached and does not interact with the flow in the bulk.
Wall-separation and vortex-breakdown zones in a solid-body rotation flow in a rotating finite-length straight circular pipe
- Zvi Rusak, Shixiao Wang
-
- Journal:
- Journal of Fluid Mechanics / Volume 759 / 25 November 2014
- Published online by Cambridge University Press:
- 24 October 2014, pp. 321-359
-
- Article
- Export citation
-
The incompressible, inviscid and axisymmetric dynamics of perturbations on a solid-body rotation flow with a uniform axial velocity in a rotating, finite-length, straight, circular pipe are studied via global analysis techniques and numerical simulations. The investigation establishes the coexistence of both axisymmetric wall-separation and vortex-breakdown zones above a critical swirl level, ${\it\omega}_{1}$. We first describe the bifurcation diagram of steady-state solutions of the flow problem as a function of the swirl ratio ${\it\omega}$. We prove that the base columnar flow is a unique steady-state solution when ${\it\omega}$ is below ${\it\omega}_{1}$. This state is asymptotically stable and a global attractor of the flow dynamics. However, when ${\it\omega}>{\it\omega}_{1}$, we reveal, in addition to the base columnar flow, the coexistence of states that describe swirling flows around either centreline stagnant breakdown zones or wall quasi-stagnant zones, where both the axial and radial velocities vanish. We demonstrate that when ${\it\omega}>{\it\omega}_{1}$, the base columnar flow is a min–max point of an energy functional that governs the problem, while the swirling flows around the quasi-stagnant and stagnant zones are global and local minimizer states and become attractors of the flow dynamics. We also find additional min–max states that are transient attractors of the flow dynamics. Numerical simulations describe the evolution of perturbations on above-critical columnar states to either the breakdown or the wall-separation states. The growth of perturbations in both cases is composed of a linear stage of the evolution, with growth rates accurately predicted by the analysis of Wang & Rusak (Phys. Fluids, vol. 8, 1996a, pp. 1007–1016), followed by a stage of saturation to either one of the separation zone states. The wall-separation states have the same chance of appearing as that of vortex-breakdown states and there is no hysteresis loop between them. This is strikingly different from the dynamics of vortices with medium or narrow vortical core size in a pipe.
On the active feedback control of a swirling flow in a finite-length pipe
- Shixiao Wang, Zvi Rusak, Steve Taylor, Rui Gong
-
- Journal:
- Journal of Fluid Mechanics / Volume 737 / 25 December 2013
- Published online by Cambridge University Press:
- 25 November 2013, pp. 280-307
-
- Article
- Export citation
-
The physical properties of a recently proposed feedback-stabilization method of a vortex flow in a finite-length straight pipe are studied for the case of a solid-body rotation flow. In the natural case, when the swirl ratio is beyond a certain critical level, linearly unstable modes appear in sequence as the swirl level is increased. Based on an asymptotic long-wave (long-pipe) approach, the global feedback control method is shown to enforce the decay in time of the perturbation’s kinetic energy and thereby quench all of the instability modes for a swirl range above the critical swirl level. The effectiveness of an extended version of this feedback flow control approach is further analysed through a detailed mode analysis of the full linear control problem for a solid-body rotation flow in a finite-length pipe that is not necessarily long. We first rigourously prove the asymptotic decay in time of all modes with real growth rates. We then compute the growth rate and shape of all modes according to the full linearized control problem for swirl levels up to 50 % above the critical level. We demonstrate that the flow is stabilized in the whole swirl range and can be even further stabilized for higher swirl levels. However, the control effectiveness is sensitive to the choice of the feedback control gain. A potentially best range of the gain is identified. An inadequate level of gain, either insufficient or excessive, could lead to a marginal control or failure of the control method at high swirl levels. The robustness of the proposed control law to stabilize both initial waves and continuous inlet flow perturbations and the elimination of the vortex breakdown process are demonstrated through numerical computations.
Stall onset on aerofoils at low to moderately high Reynolds number flows
- Wallace J. Morris II, Zvi Rusak
-
- Journal:
- Journal of Fluid Mechanics / Volume 733 / 25 October 2013
- Published online by Cambridge University Press:
- 24 September 2013, pp. 439-472
-
- Article
- Export citation
-
The inception of leading-edge stall on stationary, two-dimensional, smooth, thin aerofoils at low to moderately high chord Reynolds number flows is investigated by a reduced-order, multiscale model problem via numerical simulations. The asymptotic theory demonstrates that a subsonic flow about a thin aerofoil can be described in terms of an outer region, around most of the aerofoil’s chord, and an inner region, around the nose, that asymptotically match each other. The flow in the outer region is dominated by the classical thin aerofoil theory. Scaled (magnified) coordinates and a modified (smaller) Reynolds number $(R{e}_{M} )$ are used to correctly account for the nonlinear behaviour and extreme velocity changes in the inner region, where both the near-stagnation and high suction areas occur. It results in a model problem of a uniform, incompressible and viscous flow past a semi-infinite parabola with a far-field circulation governed by a parameter $\tilde {A} $ that is related to the aerofoil’s angle of attack, nose radius of curvature, thickness ratio, and camber. The model flow problem is solved for various values of $\tilde {A} $ through numerical simulations based on the unsteady Navier–Stokes equations. The value ${\tilde {A} }_{s} $ where a global separation zone first erupts in the nose flow, accompanied by loss of peak streamwise velocity ahead of it and change in shedding frequency behind it, is determined as a function of $R{e}_{M} $. These values indicate the stall onset on the aerofoil at various flow conditions. It is found that ${\tilde {A} }_{s} $ decreases with $R{e}_{M} $ until some limit $R{e}_{M} $ (${\sim }300$) and then increases with further increase of Reynolds number. At low values of $R{e}_{M} $ the flow is laminar and steady, even when stall occurs. The flow in this regime is dominated by the increasing effect of the adverse pressure gradient, which eventually overcomes the ability of the viscous stress to keep the boundary layer attached to the aerofoil. The change in the nature of stall at the limit $R{e}_{M} $ is attributed to the appearance of downstream travelling waves in the boundary layer that shed from the marginal separation zone and grow in size with either $\tilde {A} $ or $R{e}_{M} $. These unsteady, convective vortical structures relax the effect of the adverse pressure gradient on the viscous boundary layer to delay the onset of stall in the mean flow to higher values of ${\tilde {A} }_{s} $. Computed results show agreement with marginal separation theory at low $R{e}_{M} $ and with available experimental data at higher $R{e}_{M} $. This simplified approach provides a universal criterion to determine the stall angle of stationary thin aerofoils with a parabolic nose.
Axisymmetric swirling flow around a vortex breakdown point
- Zvi Rusak
-
- Journal:
- Journal of Fluid Mechanics / Volume 323 / 25 September 1996
- Published online by Cambridge University Press:
- 26 April 2006, pp. 79-105
-
- Article
- Export citation
-
The structure of an axisymmetric and inviscid swirling flow around a vortex breakdown point is analysed. The model assumes that a free axisymmetric bubble surface is developed in the flow with a stagnation point at its nose. The classical Squire-Long equation for the stream function ψ(x,y) (where y = r2/2) is transformed into a free boundary problem for the solution of y(x, ψ). The development of the flow is studied in three regions: the approaching flow ahead of the bubble, around the bubble nose and around the separated bubble surface. Asymptotic expansions are constructed to describe the flow ahead of and behind the stagnation point in terms of the radial distance from the vortex axis and from the bubble surface, respectively. In the intermediate region around the stagnation point, the flow is approximated by an asymptotic series of similarity terms that match the expansions in the other regions. The analysis results in two possible matching processes. Analytical expressions are given for the leading term of the intermediate expansion for each of these processes. The first solution describes a swirling flow around a constant-pressure bubble surface, over which the flow is stagnant. The second solution represents a swirling flow around a pressure-varying bubble surface, where the flow expands along the bubble nose. In both solutions, the bubble nose has a parabolic shape, and both exist only when H’ > 0 (where H’ is the derivative at the vortex centre of the total head H with the stream function ψ, and can be determined from the inlet conditions). This result is shown to be equivalent to Brown & Lopez's (1990) criterion for vortex breakdown. Good agreement is found in the region around the stagnation point between the pressure-varying bubble solution and available experimental data for axisymmetric vortex breakdown.
A transonic small-disturbance model for the propagation of weak shock waves in heterogeneous gases
- THOMAS E. GIDDINGS, ZVI RUSAK, JACOB FISH
-
- Journal:
- Journal of Fluid Mechanics / Volume 429 / 25 February 2001
- Published online by Cambridge University Press:
- 06 March 2001, pp. 255-280
-
- Article
- Export citation
-
The interaction of weak shock waves with small heterogeneities in gaseous media is studied. It is first shown that various linear theories proposed for this problem lead to pathological breakdowns or singularities in the solution near the wavefront and necessarily fail to describe this interaction. Then, a nonlinear small-disturbance model is developed. The nonlinear theory is uniformly valid and accounts for the competition between the near-sonic speed of the wavefront and the small variations of vorticity and sound speed in the heterogeneous media. This model is an extension of the transonic small-disturbance problem, with additional terms accounting for slight variations in the media. The model is used to analyse the propagation of the sonic-boom shock wave through the turbulent atmospheric boundary layer. It is found that, in this instance, the nonlinear model accounts for the observed behaviour. Various deterministic examples of interaction phenomena demonstrate good agreement with available experimental data and explain the main observed phenomena in Crow (1969).
Transonic flow of moist air around a thin airfoil with non-equilibrium and homogeneous condensation
- ZVI RUSAK, JANG-CHANG LEE
-
- Journal:
- Journal of Fluid Mechanics / Volume 403 / 25 January 2000
- Published online by Cambridge University Press:
- 25 January 2000, pp. 173-199
-
- Article
- Export citation
-
A new small-disturbance model for a steady transonic flow of moist air with non-equilibrium and homogeneous condensation around a thin airfoil is presented. The model explores the nonlinear interactions among the near-sonic speed of the flow, the small thickness ratio and angle of attack of the airfoil, and the small amount of water vapour in the air. The condensation rate is calculated according to classical nucleation and droplet growth models. The asymptotic analysis gives the similarity parameters that govern the flow problem. Also, the flow field can be described by a non-homogeneous (extended) transonic small-disturbance (TSD) equation coupled with a set of four ordinary differential equations for the calculation of the condensate (or sublimate) mass fraction. An iterative numerical scheme which combines Murman & Cole's (1971) method for the solution of the TSD equation with Simpson's integration rule for the estimation of the condensate mass production is developed. The results show good agreement with available numerical simulations using the inviscid fluid flow equations. The model is used to study the effects of humidity and of energy supply from condensation on the aerodynamic performance of airfoils.
Numerical studies of transonic BZT gas flows around thin airfoils
- CHUN-WEI WANG, ZVI RUSAK
-
- Journal:
- Journal of Fluid Mechanics / Volume 396 / 10 October 1999
- Published online by Cambridge University Press:
- 10 October 1999, pp. 109-141
-
- Article
- Export citation
-
Numerical studies of two-dimensional, transonic flows of dense gases of retrograde type, known as BZT gases, around thin airfoils are presented. The computations are guided by a recent asymptotic theory of Rusak & Wang (1997). It provides a uniformly valid solution of the flow around the entire airfoil surface which is composed of outer and inner solutions. A new transonic small-disturbance (TSD) equation solver is developed to compute the nonlinear BZT gas flow in the outer region around most of the airfoil. The flow in the inner region near the nose of the airfoil is computed by solving the problem of a sonic flow around a parabola. Numerical results of the composite solutions calculated from the asymptotic formula are compared with the solutions of the Euler equations. The comparison demonstrates that, in the leading order, the TSD solutions of BZT gas flows represent the essence of the flow character around the airfoil as computed from the Euler equations. Furthermore, guided by the asymptotic formula, the computational results demonstrate the similarity rules for transonic flows of BZT gases. There are differences between the self-similar cases that may be related to the error associated with the accuracy of the asymptotic solution. A discussion on the flow patterns around an airfoil at transonic speeds and at various upstream thermodynamic conditions is also presented. The paper provides important guidelines for future studies on this subject.